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Abstract 

We give a surface for which the Ricci Flow applied to the metric will 
increase the topological entropy of the geodesic flow. Specifically, we 
' first adapt the Melnikov method to apply to a Ricci Flow perturbation 



and then we construct a surface which is closely related to a surface of 
revolution, but does not quite have rotational symmetry. This is done 
by adapting the Liouville metric representation of a surface of revolution. 
The final steps of the Melnikov method require numerical integration. 



' 1 Introduction 

' Let M be a closed surface with cotangent bundle r : T*M M. Let oj be the 

. canonical symplectic structure and let the Hamiltonian be 

\0 ■ 1 

O ■ H{x,p) := -gx{p,p), 

. where g is the dual of the Riemannian metric. The flow 0t associated to the 

^ ' Hamiltonian system {T* M,uj, H) is called the geodesic flow of {M,g) TT. 

^ . An important dynamical invariant for a geodesic flow is the topological en- 

tropy a measure of how rich the dynamics within the flow is. On a surface 
the entropy is non-zero if and only if there exists a transversal homoclinic point, 
' and in this case there will be exponential growth in the lengths of periodic orbits 

cni. 

The Ricci Flow of a metric g is a smooth family of metrics gt satisfying 



dgt 
dt 



— = -2RiCgj, with .go = .9- 



Introduced by R. Hamilton |6j, the Ricci Flow is in some sense an analogue of 
the heat equation for the Riemannian metric [HI p. 2]. Loosely speaking, the 
Ricci Flow acts to diffuse the curvature across the manifold - homogenising the 
Riemannian structure. 

The Ricci Flow is particularly well understood on surfaces, where it simplifies 
the metric towards a constant scalar metric in a controlled fashion, [Sj |S] . 
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In JT] , Manning used a result of Katok, Knieper and Weiss PP on Anosov 
flows to show that as the metric on a surface of everywhere negative curvature 
evolves under the Ricci Flow, the topological entropy of the geodesic flow is 
strictly decreasing. 

Decreasing entropy corresponds to a simpler dynamical structure, and at a 
first instance this is what one might always expect as the Riemmanian struc- 
ture is being simplified. However, here we give an example of a surface for 
whom the Ricci Flow will increase the entropy: 'simplifying' the metric can 
break symmetries in the overlying symplectic structure and induce complicated 
dynamics. 

For ease throughout, we call a Riemannian 2-manifold a surface, and a com- 
pact surface whose geodesic flow is integrable and has a homoclinic connection 
a good surface. We let the curvature of our good surface be K, let (pt be the 
geodesic flow, let {z^, a:^,pi,p2} be standard local co-ordinates of the cotangent 
bundle, let F be the homoclinic connection, and let F be the second integral of 
motion. 

When a good surface is perturbed through metrics that satisfy the Ricci Flow 
equation, one can tell whether or not the geodesic flow still has zero entropy by 
applying the following theorem: 

Proposition A. We define 



If M{z) changes sign then the Ricci Flow will strictly increase the entropy of 
the flow. 

The Normalised Ricci Flow is a modification of the Ricci Flow that preserves 
the volume of the manifold; 



where A is the average curvature and n — dim(Af). If, in proposition^ we 
were to use the Normalised Ricci Flow instead of the Ricci Flow the conditions 
for increasing entropy would be the same. This is unsurprising as the Nor- 
malised Ricci Flow is a space-time reparametrization of the Ricci Flow, and the 
perturbative Melnikov method is the backbone of the proof. 

Whilst many examples of good surfaces were known classically - such as 
ellipsoids and poisson spheres - explicit calculation of the integral in these cases 
has proved intractable for us. We show how to construct new examples of good 
surfaces for which the hypotheses of proposition ^ hold, and go through the 
calculations in a specific case: 




Theorem B. If a torus has a Liouville metric 



where g{y) is a smooth 27r -periodic function strictly greater than —1/2 and equal 
to 2/^(1 — y^) on [—1, 1], then the topological entropy of its geodesic flow will be 
increased by the Ricci Flow. 

Since the surface above is a torus, the resuhs of Hamihon |H] imply that as 
time under the Ricci Flow tends to infinity the metric tends to the flat metric. 
Interestingly then, this is a surface for which the metric has zero entropy for 
t = Q and t — > oo, but is not always zero. 

Using a branched cover construction, we outline a related example with genus 
zero. 

Corollary C. There exists a sphere whose geodesic flow has zero entropy ini- 
tially, but the metric is such that under the Ricci Flow this entropy will increase. 

In subsection 12 . II we outline the Melnikov method and then, in 12. 21 apply it 
to a Ricci Flow perturbation to prove proposition 1X1 

Section |31 summarises the argument for constructing a surface on which 
proposition fXl gives a positive result. The properties of Liouville surfaces are 
discussed in section 13.11 as these generalise surfaces of revolution. Surfaces of 
revolution can easily be given a homoclinic connection but fail to be useful for 
another reason, as outlined in remark [3. 71 With a well chosen modification we 
construct a suitable surface in 13. 21 

The last steps involve a numerical integration and the standard graphs show- 
ing convergence are given. 

2 Melnikov's method and the Ricci Flow 
2.1 The Melnikov method 

Assume the 4-dimensional Hamiltonian system (M, oj, Hq) has a heteroclinic or 
homoclinic orbit. We will consider the effect of perturbing the Hamiltonian by 
studying the dynamics in the family of conservative systems {M,oj,H^), where 

H, = Ho + eHi+0{e'). 

We assume that the unperturbed system is integrable, with second integral 
F. We pick a homoclinic orbit and call it F. 

Remark 2.1. All arguments also apply in the case of a heteroclinic rather than 
a homoclinic connection, but for ease of exposition and without significant loss 
of generality we assume all connections are homoclinic. 

Melnikov's idea was to measure the transversal distance between the stable 
and unstable manifolds in the perturbed system, using F as a reference point. 
If the stable and unstable manifolds do not coincide then we use theorem 5 of 
Robinson for analytic Hamiltonians to show that the topological entropy 
of the system is positive for all e > sufficiently small. Alternatively, if we do 
not want to insist in our system being analytic in the relevant region we could 
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also use the results by Burns and Weiss 0]. In this exposition the dynamics in 
the regions on which we apply the Melnikov method are analytic. This method 
is very well understood: see Robinson ^S] and Wiggins which include 
discussions on why the integrals below converge. 

Let Y be the first order change in the Hamiltonian vector field. Precisely, Y 
is such that if we let Xz be the vector field associated to a Hamiltonian Z then 

Xh, = Xh, +fY + 0{e^). 

Definition 2.2. The Melnikov integral is 

/oo 
[YF)[c^t{z))dt. 
-oo 

Theorem 2.3 (Melnikov |12|. Robinson |15p. // the Melnikov integral 
M{z) changes sign, then for all e > sufficiently small the entropy of (Af , lu, H^) 
is positive. 



2.2 Proof of proposition 1X1 

We let {M,gt) be a Riemannian two-manifold solution of the Ricci Flow. Since 
the Gaussian curvature satisfies K = 2Ric (also a function of t), 

We refer to the standard literature for results: in particular that the Ricci 
Flow of a surface exists for t S [0,T], with T > 0. For more information, 
including the recent proof of the Geometrization Conjecture, see Chow and 
Knopf inj and Morgan and Tian 

We consider the family of geodesic flows of (M, gt) given by the Hamiltonian 
systems {T* M,uj, gtip,p)) for t e [0,r], T > 0. 

Assume the geodesic flow at t = is integrable with second integral F and 
that the flow has a homoclinic orbit. We look to apply Melnikov's method to 
see whether the initial Ricci Flow breaks the homoclinic orbit and hence the 
perturbed system has positive entropy. 

Because the metric is being flowed on the tangent bundle but we need the 
norm of covectors the following lemma is necessary. 

Lemma 2.4. 

dgt ^ . 

Proof. We know gtgt is the identity tensor. On differentiating and applying the 
product rule the result is immediate. □ 

Remark 2.5. The Hamiltonian is a function with domain T*M, and as we 
will see the Gaussian curvature arises naturally when we consider first order 
perturbations. However, the Gaussian curvature has domain M; therefore we 
need to use K :— K o t, where r is the footpoint projection mentioned in the 
introduction. 
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Lemma 2.6. For a Ricci Flow perturbation of the geodesic flow, Y = Xj^^^. 
Proof. 

1 1 € d 

z 2 I at t=o 

= Ha{p) +eKHo{p) +0{e^). 

Therefore the perturbed Hamiltonian vector field takes the form = X^o + 
^-^KHo + '^(*^^)' which imphes that Y = X^^^. □ 

We use this to simphfy YF in theorem 12.31 Using the poisson brackets 
induced by to 

YF = X^^F = {F,kH^} 

= K{F,Ho} + Ho{F,K} 
= Ho{F,k} 

since F is an integral. 

Proposition 2.7. In local co-ordinates {x^ ,pi,p2} on T*M, 

dK dF 
YF = . 

dx^ dpj 

Proof. Since we can consider the flow on the level set Hq ^ 1 
Ho{F,k} = {F,k} 

= -dkiXp) 

= -dKodTiXp) 

= dK{cp), 

where (p = —dT{Xp). The proposition will be proved if we can find ip in local 
co-ordinates. We have 

dF dF 

but as these co-ordinates are Darboux, 

^ _ dF^d_ _ dF d 
^ dx^ dpi dpj dxi 

The result follows on applying dr, i.e. by setting the components to 
zero. □ 

Together with the result of Robinson ^S], on how the Melnikov integral 
implies positive entropy, we have proved proposition 1X1 

Proposition A. We define 

f°° dK dF ^ , 

If M{z) changes sign then the Ricci Flow will strictly increase the entropy of 
the flow. 
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3 Construction of a tractable good surface 



Definition 3.1. A Liouville surface is a surface whose metric is locally of the 
form {f{x) + g(y))(dx'^ + dy^) for some /, g real functions with f + g positive. 

All the good surfaces mentioned in the introduction were quadratically inte- 
grable. It is a classical result (see Darboux [7]) that such surfaces are Liouville 
surfaces. The example surface of this paper is also a Liouville surface. 



3.1 Simplifications on a Liouville surface 

Throughout we let A — f{x) + so that ds^ = A{dx^ + dy^). Our standard 
reference is A. V. Bolsinov and A. T. Fomenko, |2|. 

The Liouville metric is an isothermal parametrization. Hence the Gaussian 
curvature is 

/.^-A-AlogVA.l^-^, (1) 

and the partial derivatives are 

= {-f"'A^ + 2f'g"A + Aff'A - 3f'g'^ ~ 3/'3) /{2A'^), 
Ky = {-g"'A^ + 2f"g'A + Ag'g"A - if'^g' - ig'^) /{2A% 



(2) 



Lemma 3.2. The geodesic flow on a Liouville surface is integrable, with integral 
F=pl-fH. 

Proof. By direct calculation. □ 

Lemma 3.3. For any geodesic 7 = (x, y) on a Liouville surface, there exists a 
constant a such that 

(3) 

Proof. We pass to the tangent bundle, using = Ax, Py = Ay. Since H and 
F are preserved by the flow, so is 

F 9pI - fpl gx'' - fy' 




Now we solve for y/x. □ 

To simplify the exposition we throw away some generality at this stage. With 
theorem IbI in mind, we assume the Liouville metric is defined on a cylindrical 
neighbourhood arising from a lifted co-ordinate strip R x [—a, a]. 

A geodesic on the surface lifts to a line in the co-ordinate strip. Hence a point 
z{t) moving along a geodesic on the surface determines a point {x{t),y{t)) e 
R X [—a, a]. Once we set one point in the co-ordinate chart, x{t) and y{t) are 
uniquely defined. If we assume y/x is bounded, then each value of x uniquely 
determines a value of y, and we can equivalently describe the geodesic by the 
function y ~ y{x). 
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We make sure that geodesies under consideration are trapped in this neigh- 
bourhood, and that the derivative bounded, in remark [3 . 81 below. 

Corollary 3.4. The Melnikov integral in proposition simplifies on such a 
Liouville strip; 

r°° dK OF dv 
J ^Q^Q-l(M^))dt = 2j ^K^{x,y)g{y)-Ky{x,y)f{x)^dx. (4) 

Proof. We substitute F — p^ ~ fH as a second integral. 
f°° dK dF ^ ^ , , r°° ^K^g{y)Px~ Kyf(x)py , ^ , , 

/OO nOO 
K^g{y)x{(pt{z))dt-2 / Kyf{x)y{(pt{z))dt. 
-oo <y —oo 

Then we make the substitution x = x(t) and view y as a function of x as we 
run along the orbit (t>t{z). The improper integrals must all converge absolutely 
for the argument to hold, but this has been assumed since the beginning. □ 

For clarity we define I{x,y) Kx{x,y)g{y) - Ky{x,y)f{x)dy/dx. 

Remark 3.5. From now on we restrict the flow to H^^{1), a 3-manifold. Ob- 
viously local co-ordinates are {x, y, 6}, where tan^ = p^/py 

Proposition 3.6. In addition to our current hypotheses, if g'{0) = and 

ff"(0)>r(0)//(0)-/"(0) 
then y = is a hyperbolic geodesic. 

Proof. Liouville's theorem implies the phase space will foliate into 2-manifolds 
wherever dF and dH are linearly independent. Direct calculation shows that 
they are nowhere parallel; the foliation fails only because dF = 0. 

This foliation cannot hold true on y = if it is a hyperbolic geodesic, and 
so the proof will follow on considering conditions for dF = 0, which are 

-/' sin^ edx + g' cos^ 9dy - (/ + g) sin 29de = 

The last coefficient, (/ + g) sin 20, is only zero when sin 20 is zero, a.s f + g 
is always positive. Without loss of generality this implies that 6 = tt/2. The 
second term yields g'{Q) — 0. 

As 9 = tt/2 at every point along y = 0, the axis is a geodesic. By substi- 
tuting the condition for ^"(O) into ^ we see that the curvature along y = is 
everywhere negative. Hyperbolicity of the geodesic along the axis then follows 
from a standard Jacobi field argument. □ 
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Remark 3.7. It is well known that the Ricci Flow preserves the isometry group 
of the metric, see Hamilton 6, p. 3]. Thus if we were to Ricci Flow a surface of 
revolution, the symmetry would be preserved and the surface would remain a 
surface of revolution, and so remain integrable with second integral the Clairaut 
integral. 

This implies the geodesic flow on a surface of revolution will remain inte- 
grable along the Ricci Flow; on surface flows this implies that the entropy re- 
mains zero. Hence if there was a homoclinic connection in the geodesic flow, no 
homoclinic tangle would form: the Melnikov integral on a surface of revolution 
should always be zero. 

A surface of revolution is exactly a Liouville surface with / = 0. By the 
above argument, for a surface of revolution the hypotheses of proposition ^ 
should not apply and the Melnikov integral of corollary 13.41 should always be 
zero: 



M{z) = 2 J K^{x,y)g{y)~Ky{x,y)f{x)^dx 
= 2 J 0.g{y)-Kyix,y)-0-^dx = 0, 

as expected. 

In the next subsection we look to break the symmetry by modifying 

3.2 Explicit calculation of the Melnikov integral 

On a good surface, we have a hyperbolic orbit c with a homoclinic connection 
F. The geodesic flow is integrable with second integral F. 

To calculate an explicit local relation y = 7(2;) which describes a geodesic in 
F we must solve the orbit equation in lemma By continuity, the value of a 
on the homoclinic connection is the same as on c, the hyperbolic orbit. We will 
assume that on c we have y = and Pj, = 0, as in proposition Since we are 
on the unit tangent bundle these assumptions imply that a = g(0). 

We set g{0) = and from © obtain the following equation: 



The sign determines whether we are considering the trajectory above or below 
2/ = 0; the geodesies y = ±7(0:). Due to the symmetry of the example outlined 
below we consider the positive solution without loss of generality. 

We look to model our surface on a surface of revolution, as in remark l3. 71 
and wish to cover a cylindrical open set of a surface with Liouville co-ordinates. 
We will need / periodic in order to take a quotient of R x [—1,1] and form a 
cylinder (see figure P). 
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hyperbolic orbit 



stable orbits 



Figure 1: Parametrization of the surface by a horizontal strip of M x [—1,1]. 
The local co-ordinates must be periodic in x in order to be able to quotient to 
a cylinder. 



To ensure / is always positive after periodic modification, we take our Liou- 
ville surface with / = 1. 

Remark 3.8. If / is constant, then the metric, function 7, curvature, and g 
are all dependent: in particular, we can set 7 rather than set g. Hence we can 
ensure that the geodesic is restricted to the co-ordinate chart, M x [—1, 1], and 
that a value of x uniquely determines a value of y along a geodesic in T . 

We use this observation to construct the surface, and then perturb / away 
from 1. Although various functions are suitable candidates for 7, we chose 
7(2;) = 1/2 sech X as the convergences are exponential and the surface is smooth. 
This induces g{y) = ?/^(l — y^). As expected, this is a function with a single 
minimum - exactly as in the surface of revolution case, remark [3.71 

We choose our perturbation of / by asking what properties we require in 
y — ^{x). From the form of g determined by 7 in remark f3. 81 and from the 
orbit equation, 



for some constant of integration k. The function / must be smooth and periodic 
- and hence bounded. The simplest candidate is f{x) — (H-/icosa;)~^, for some 
small /i, so that a geodesic in the homoclinic connection F is determined by 



Remark 3.9. The Melnikov function is evaluated sX z = {x,y), a point on the 
homoclinic orbit. Since k fixes z and vice versa we consider M{k) instead of 
M{z). In the particular example just outhned M{k) is a an odd function of k: 
the aim is simply to find a value of k such that M{k) is non-zero. 




y — 7(k + x + iisinx) = — sech (k, + x + /isinx). 



(6) 
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Unfortunately, attempts to show that integrals of this form are non-zero 
through analytic means alone have so far been unsuccessful. We proceed with 
a method that requires numerical analysis for one part of the argument. 

We split the integral into two parts, M(k) = Al{k) + Bl{k), where 

Al{k) — / I {x, ■f {k + X + fj, sin x))dx 

and Bl(k) is the integral over the remaining outlying regions. 

The argument proceeds as follows. We fix k and L. We show that Al is large 
to within some accuracy e, whereas \Bl\ is small. Then AI{k) > Al — b—Bl > 0, 
so that M{k) is non-zero and we are done. 

Lemma 3.10. For < /i < 1/4, 

Bl{k) < 1200(1 + 2e2l«l)7(L)2. 

Proof. 

/(x, y) = K^g - Kyfdy/dx = fdy/dx(K^dy/dx - Ky). 

f is bounded by (1 — /i)^^ and dy/dx goes to zero exponentially, so it remains 
to bound the change in curvature. 

As an example we bound Kx along 7; Ky follows in a very similar manner. 
We bound /, g and their derivatives before plugging them into Q for and 
simplifying. 

\f{^)\ < (1 ^ mentioned, and |/'(a;)| < 2/i(l — n)^^ is also trivial. 

„ 2/x cos^ a; — cos a; — 3/i 

^ " (l + Aicosx)4 • 

Given < fj, < 1/4, it is simple to show that the norm of the numerator 
has maximum 1 + fi. This gives |/"(a;)| < 2^(1 -I- ^)(1 — Analogously 

<2M1 + 4a*)(1-m)-'- 
Since < /i < 1/4, we have 

\f{x)\<2, |/'(x)|<2, 
|/"(a:)|<2, \f"'{x)\<5. 

Consider the half x > 0. We know that here y is monotonically decreasing 
on X. The derivatives of g{y) — j/^(l — y^) can be bounded by plugging in the 
maximum value of y, which is less than ^{L — \k\ — fi) = 70. This 70 is less 
than 1, as j{x) < 1 for all x, and we use this below. The argument for x < is 
symmetric. 

Thus 

l5(y)l<7o^ < 270, 

\g"{y)\<2, |.9"'(y)|< 247o. 

We note A'^ = (/(x) + giy))~^ < (1 + A^)^- On substitution into 

\Kx\ < 2470(1 + (4 + 8)(1 + m)*' + (670' +47o')(l + A^)' 
< (1 + Ai)^(24 + 12 + 10) < 276. 
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The calculation of \Ky\ < 28870, which tends to zero as |i + K+ l| 00, 
reflects Ky{x,0) = 0. 

Lastly, \dy/dx\ < 5j{L)/A as tanh is bounded by 1, and another easy calcu- 
lation gives 7o < '^{L) ■ exp(2/i + 2|kP). 

We are now in a position to show 



I{x, y)dx 



< 



1/(^)1 \K: 



dy_ 

dx 



dy 
dx 



dx 



< 2(2767(L) + 2887o) 

< 600(1 + 2e2|'"l)7(L)^ 



7(00) 



-dy 



7(i) 



Bl is the sum of two integrals both bounded by this value, which completes 
the analysis. □ 

Remark 3.11. This bomid is very rough - perhaps to the casual observer 
recklessly so - but it won't matter as 7(2:) — > exponentially fast. For example, 
with |k| < 3 and < ^ < 1/4, we only need L = 10 for Bio{k) < 0.005. 

We use Maple to implement the fourth order classical Runge-Kutta method 
on Al{k), with k = 1 and (arbitrarily) /i = 1/8 7^ 0. Evidence that the error 
is being controlled as usual is shown in figure [21 which shows convergence as 
expected. 




0.5 1 1.5 2 2.5 3 

ln( Step Size ) 



Figure 2: loge-loge graph of convergence as step size decreases, with k — I and 
/i — 1/8. Standard convergence is given as a reference gradient. Initial noise is 
due to inaccurately low step sizes. 
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-1 -0.5 0.5 1 1.5 2 

Figure 3: Melnikov Integral as a function of k (a;-axis in multiples of tt), with 



Table 1: Values of ^l(I) and the bound for Bl(1) against L, to 8 significant 
figures, with ^ = 1/8. 



L 


Al{1) 


bound for Bl{1) 


1 


1.5990894 


7.951683E4 


2 


1.0073283 


1.337684E4 


3 


0.59562584 


1.868009E1 


5 


0.48440447 


3.438049 


10 


0.47929061 


1.561G13E-4 


20 


0.47929047 


3.217488E-13 



Remark 3.12. As expected, if we hold L constant and calculate Al{k) in terms 
of varying k (with a reasonable accuracy) we get a sinusoidal graph, figure 01 

Tabled shows how the values of Al{1) converge as L increases, and we also 
include the bound for Bl{1) calculated in lemma mill This gives irrefutable 
evidence that the Melnikov integral is non-zero: 0.479 to three significant figures. 
By proposition we have shown 

Theorem B. If a torus has a Liouville metric 

where g{y) is a smooth 2'K-periodic function strictly greater than —1/2 and equal 
to j/^(l — y^) on [—1, 1], then the topological entropy of its geodesic flow will be 
increased by the Ricci Flow. 
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Proof of corollary C. Quotienting by the involution a{x,y) — {—x^—y) gives 
a branched double cover of over S^. We assume F and G are both even 
functions, so they are invariant under a. If F and G are chosen such that 
F(x) + G{y) vanishes at the four branch points in a suitable manner, then the 
induced symmetric tensor on S"^ will be a Liouville metric. 
The 'suitable manner' is that near the zeroes 



F{t) + G{it) =0, 



(7) 



and that at any zero F has a non-zero second derivative. These relations arise 
since the branch points are double branch points, as in 0. 

We take F similar to / above, but subtracting a constant so that it has 
zeroes; F{x) = (1 + iicosx)~^ — (1 + By adding the same constant to 

G the metric is unaffected and the argument for the Ricci Flow to increase the 
entropy still holds. We shift the strip on which the argument was applied away 
from the branch points using y — y — Z. Hence a suitable function for G is a 
smooth function reflecting at and 3 such that 



G{y) = {l + ^l)-^ + 



See figure 0] 



yHl-y^) ye [2,3] 
-(l + ^coshy)-2 ye [0,1] 



Toms in the plane 



(reflected) surface of l evolution-like region 



Branch points 




cc shine-like region 



Dynamics as required in mis strip 



G(y) 



Figure 4: How F and G are adapted to the branched cover of over . Zeroes 
of F and G occur simultaneously at the branch points. 
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Thus wc have defined a Riemannian metric on ]R^/(27rZ® 12Z)/(T, which is 
topologically a sphere. On this surface the geodesic flow has zero topological 
entropy, but under the Ricci Flow this entropy will become strictly positive. □ 
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